Calculating model parameters 

This supplementary material presents an algorithm for calculating the param- 
eters of the network model for a circuit of resistors. The network equations 
describing such a circuit with N nodes connected to a dc power supply are, 



for i,j = 1,2, ... ,N. Here, A is an NxN matrix, and V^ is an TV-element 
vector containing the node potentials of the unperturbed circuit. The defin- 
ing feature of the network model is that these constants are determined by 
experiment. 

In order to calculate the parameters, N + 1 experiments need to be per- 
formed. Measuring the node potentials for the unperturbed circuit directly 
determines V^ . The remaining experiments involve fixing the potential of a 
single node and measuring the potentials of the TV — 1 remaining nodes. This 
is done for each node. In this work, we set the fixed potential to be zero; 
however, this is not required, and any convenient potential can be used. 

The procedure for calculating the model parameters will be demonstrated 
on a circuit with three nodes for simplicity. It can be extended easily to larger 
circuits. During a measurement of the unperturbed circuit, the equations 
describing the node potential values are simply Eqs. (1), i.e., 

= A u (Vx - V; (o) ) + A 12 (v 2 - K, (o) ) + A 13 (v 3 - vf >) (2a) 

o = a 21 (Vx - v; (o) ) + a 22 (v 2 - \/ 2 (0) ) + a 23 (y 3 - v 3 (o) ) (2b) 

= A 31 (y x - Vi 0) ) + A 32 (y 2 - V 2 (o) ) + A 33 (V 3 - \/ 3 (o) ) . (2c) 

Clearly, the potentials V^ satisfy these equations. 

The remaining experiments involve fixing a single node potential to zero 
by grounding. When node one is grounded, for example, the equation de- 
scribing its behavior is no longer valid because the node is under external 
control and no longer influenced by the others. In such a situation, the node 
potentials of the remaining nodes will be modified from their normal val- 
ues. The new potentials, when measured, are V 1 ^ 



t/(1) yd) t/(1) 
v l ■> v 2 ■> v 3 



with 



Vi = 0. The equations corresponding to such an experiment can be written 
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Next, node two is grounded to obtain the potentials V^ 2 ) and the equa- 
tions 



= A n (v; (2) - V; (o) ) + A 12 (V 2 (2) - V 2 ^) + A 13 ( 
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Grounding node three generates, 
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These are all the equations needed to calculate the model parameters. We 
can now collect the equations containing similar elements of A. For example, 
Eqs. (2a), (4a), (5a), all contain the first-row elements. These can be written 
in matrix form, i.e., 
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The matrix contains the data collected from experiment. The unperturbed 
data has been included in the first equation for the moment. If a computer, 
e.g., with software such as Matlab, is used to solve this system of homoge- 
neous equations for [An, A 12 , A 13 ], it will return a vector of zeros, which is 
not the desired solution, of course. A few manipulations are necessary to 
avoid this. First, the V{° terms in the left column are moved to the other 
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side of the equations. The remaining elements of the top row are simplified 
by inserting a value of zero. This results in, 
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The first equation produces A u = 1, and the remaining parameters can now 
be calculated. 

A similar procedure can be followed to calculate the remaining rows of A. 
As discussed in the main text, for a circuit of resistors, the equations gener- 
ated by the network method are the same as those obtained using Kirchhoff 's 
laws. 

We have included, in the Supplementary Materials, a Matlab function 
which automates the calculation of the A matrix. It requires a two inputs, the 
total number of nodes (N) and the matrix of data obtained by experiment. 
The matrix should have the form, 
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Note the position of the unperturbed data. It is included at the bottom of 
the matrix. 



